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In this note we employ concurrently techniques of generalized zeto-functions and compactification 
methods introduced in previous publications, to study the Halperin-Lubensky-Ma theory of induced 
weak first-order phase transitions applied to type-I superconducting films. We obtain closed formulas 
to the critical temperature and to the size temperature as functions of the film thickness. 
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The Halperin-Lubensky-Ma (HLM) effect appeared 
about three decades ago [S| . It predicts a weak first-order 
phase transition in superconductors. This fact emerges 
by considering in the Ginzburg-Landau (GL) model the 
interaction between the intrinsic magnetic fluctuations 
and the order parameter. More generally, all the phys- 
ical systems that are described by the Abelian-Higgs 
model present this phenomenon. Some examples are: 
the nematic-smectic A phase transition in liquid crystals 
0, 0, with the interaction between the smectic scalar 
order parameter and the vector director; the massless 
scalar electrodynamics in field theory .3], which presents 
the gauge field acquiring mass as an effect of its coupling 
with the scalar field. 

However, the size temperature associated to the HLM 
effect was determined to be too small which makes it 
very difficult to be detected experimentally. Recently, 
the mentioned effect has been studied in the context of 
type-I superconducting films 0, IE IE S Hi ■ In particu- 
lar, Refs. p], IE El suggest the enhancement of first-order 
transition in superconducting films with respect to that 
in bulk materials, which at least qualitatively is corrob- 
orated in Ref. Q- 

In order to have a better understanding of HLM effect 
in films, in this note we extend some questions raised 
in Ref. In a field theoretical approach, we con- 

sider the GL model submitted to confinement between 
two planes a distance L apart from one another. This 
is done using a spatial compactification formalism pre- 
sented in recent works [E El El 0, El- Physically, 
for dimension d = 3, L corresponds to the thickness of a 
film-like superconducting material. We take into account 
the gauge fluctuations in absence of external magnetic 
field, and in the approximation of uniform order param- 
eter ip (x) = const. We investigate the critical behavior 
of the system as a function of the film thickness L, and 
in particular we focus on the L-dependence of the inter- 
esting thermodynamical quantities, as well as we discuss 
the plausibility of the presented results. 

Let us consider the Hamiltonian density of the GL 
model in Euclidean rf-dimensional space, 



-™olVf 



(i) 



where ip = ip (x) is a complex field, and ?71q is the 
bare mass and = (x) {u,v = l...d) is the gauge 
field. Notice that we are working in the mean field con- 
vention |14| . Accordingly, we use natural units (h — 
c = I), and employ £o (the intrinsic coherence length) 
and KbTq (Tq corresponds to the bulk mean field tran- 
sition temperature), as length and energy scales, re- 
spectively. The fields and coordinates are rescaled by 
ip -* ipnew = v^oAs^b^, A -> A new = y/£a/k B T A 
and x — > x new = x/^q. Thus, the Hamiltonian density 
in Eq.Q is dimensionlcss. Also the coupling constants e 
and A are dimensionless, as well as the bare mass, which 
is given by m 2 , = T/Tq — 1. 

The definition of the field as ip new — 4> e , together 
with the gauge transformation A new — » a— 1/eVfl, allows 
us to work with Ti. manifestly gauge invariant [l4l Hi} . 
Then, the integrations over ip new , tp^ ew and A new in the 
generating functional change to integrals over <f>, a and 
9. Since the functional integration over 9 is Gaussian, 
it can be performed directly. Taking explicitly the ap- 
proximation in which the scalar field is spatially uniform, 
that is (j) ( x ) ~ 4> = const., the functional integral over 
a can be done, yielding the following free energy density 



T= l -ml^ + ^ + V{ ( j ) ), 



(2) 



where 



H 



1 



(d^-d^f + Kd^-ieA^f 



v{<j>) = — in { [(-v 2 + c 2 4> 2 ) <v + 5 ( x - y)} ■ 

is the contribution coming from the a-integration. In the 
case that the model has no restrictions of spatial con- 
finement, it can be shown that performing the trace in 
the right hand side of Eq. 10 and using the momentum 
space representation, a 3 -dependent term is generated 
for d — 3, which suggests that systems described by the 
free energy @ undergoes a first-order phase transition 

Q 

However, for situations in which the model is restricted 
to boundary conditions, Eq. @ must be treated with 
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some additional techniques. Let us consider the system 
confined between two parallel planes, normal to the x^- 
axis, a distance L apart from one another (we will use 
the dimensionless form I = L/£q). We use coordinates 
x = (z,Xd), where z is a (d — 1) -dimensional vector, 
with corresponding momenta k = (q, kd). Accordingly 
to 0, in this confined situation, the term corresponding 
to V(4>) in Eq. J3J should be written as 



v(<M) 



r/(0; <M), 



(4) 



where the prime means derivation with respect to the 
first argument in r\. The function r\ is associated with 
the eigenvalues of the operator given by Eq. with the 
compactification of one dimension, 
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We can perform the integration in the above equation 
with the help of dimensional regularization techniques; 
we get 



T)(s;c,l) = 



T(s) 



d-1 
s — ; wi 



(6) 



where we have defined 



w N n N 



{n; } — — OO 

(7) 

with Wi = h/h (i = 1,...,N) and c 2 = (ecj)l/2Tr) . 
N stands for the number of compactified dimensions, 
in the general case of a A^-dimensional compactified 
subspace(./V < d). For us it is TV = 1, and we identify 
h = I. 

As is shown in Ref. |9(, Eq. (JSJ and consequently Eq. 
(0J can be analytically extended to all values of d. For 
odd ^(remember that we are interested in d = 3 with one 
compactified dimension), the analytical continuation of 
V((f>,L) has the following expression[!| , 



(8) 



where p = For small values of c 2 , c 2 << 1, it is 

2 

possible to use the binomial expansion for Zf in order 
to expand it in powers of the field 0, 



Zf (q; wi) = c- 2 * + Y, ^-Tiq, j)Z x (q + j; Wl ) 



3=0 



where 



r(g + j) 
r(«) : 



(9) 
(10) 



and 



Z N (v,{wi})= Y ' [wlnl + ... + w 2 N n%] ".(11) 

{n j } — — oo 

The prime in l|llf) means that the term {rii} — is ex- 
cluded from the summation. For N = 1, w\ = 1, and it is 
easily seen that Z\ (y,W\ = 1) is related to the Riemann 
zeta-function, C,( v ) by, 



Zi (i/,wi)=2C(2i/), 



(12) 



Then the derivative of ZJ 2 (g; t«i = 1) for c 2 << 1 in Eq. 
© above with respect to q gives the following expression, 



Z' x c2<<1 {q; ioi = 1) = -2c- 29 In c + 2 ]T c ; 



2j 



3=0 



x {T'(g, j)C (2? + 2j) + 2T(g,i)C' (2<? + 2j)} . (13) 

We take into account in the expansion in Eq. (|13fl only 
the c 2 -dependent terms up to the second order, which is 
consistent with the spirit of the GL model. The terms (j> 6 , 
<fi 8 , are considered as irrelevant in the neighborhood 
of the transition. So, Eq. JSJ becomes 



V(c,l) 



-^y£r{- 2c2pinc - 2c2 [C(^ + 2) 

2K'(-2p + 2)] + c 4 [(-2p + l)C(-2p + 4) 
2p(-p + l)C'(-2p + 4)]}. (14) 



Hence, considering the particular case of our interest, 
d = 3 [p = 1), and writing c explicitly in terms of 4> a- n d 
i, we obtain the correct expression for the effective free 
energy density, 

^;O^mV + ^ 4 -W> 2 ln0, (15) 

with the following definitions of the coefficients m 2 , u 
and v, 



m 2 = m 2 + v (1 — 2 IneZ) , 
e 4 Z 



u = A 



247r' 
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(16) 
(17) 

(18) 



Notice that the concavity condition for the free energy 
(|15|l restricts the model to strictly positive values of u . 
This implies, using Eq. I|17fl . in the upper limit for the 
film thickness, 



jmax 



24ttA 



(19) 



It is also worthy mentioning that Eq. i|15fl above has the 
same formal dependence on the order parameter <f> as in 
Ref. Eq. I)15|l describes a first-order transition in 

systems in a form of a film. However, we have obtained 
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a logarithmic /-dependence in Eq. (|16f) above for the 
coefficient of the </> 2 term, as well as a linear dependence 
on / for the self-coupling constant u, given by Eq. l(T7|. 

In order to have a better understanding of the /- 
dependence of this phase transition, we can study the 
behaviour of the system from Eq. Ijl5(l , in a similar way 
as is done in Refs. and Q. Taking the first deriva- 
tive of J- with respect to 0, we have T' = 4>g, where the 
function g = g(4>) is 



fJ 



v + — (j) — v<p In 0. 



(20) 



Therefore, for the equation of state T' = 0, we get the 
following possible solutions: <fi = </>oa' = 0, which de- 
scribes the normal phase, and positive solutions coming 
from g = 0, which correspond the broken phase. We 
denote them by 4>qb- Furthermore, the solutions must 
produce minima in J- to be stable, i.e. they must obey 
the condition !F"(<po) > 0. Then, since T" = g'cj) + g, 
we clearly see the normal phase is always stable. On 
the other hand, in the superconducting phase we have 

4>ob > and g'(<poB) > 0, which yields O w > J~^- 
Notice that these solutions are real and finite for u > 
(which as already remarked is also the condition of con- 
cavity of the model) . In addition, we can obtain a par- 
ticular value of 4>qb for which the free energy density 
vanishes. Writing T in terms of g, given by Eq. I|2(J|) . we 
obtain for J-(4>qb) = the value c/)qb = 4>ob,e — ^\[\- 
Notice that there are different values, 4>qn aud (/)ob, E , 
which imply in T = 0; this situation corresponds to the 
equilibrium transition point of the first order transition, 
as remarked in Ref. 

Thus, it is possible to obtain the value of the coefficient 



at the equilibrium point, J-(4>ob,e) — 0; we get 



v I In ■ 



-it; 



1 



(21) 



So, from Eq. I|16|) an d expliciting the dependence of u 
and v in /(= £,o/L) we get the following expression for 
the /-dependent equilibrium transition temperature, 



T E (l) = T 



1 



4ttZ 



In 
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(22) 



To obtain an quantitative estimate from Eq. I)22fl . let us 
remember the 3-dimensional expressions for the coupling 
constants Il4l. 



A « 111.08 



T F 



2.59 



(23) 



where Tp, vf and a are respectively the Fermi tempera- 
ture, Fermi velocity and fine structure constant (remem- 
ber also / = £ /£)). 

Considering as example an aluminum sample, and tak- 
ing their tabulated values, T = 1.19°^, T F = 13.6 x 
10 4 K, vp = 2 x 10 6 to/s and £ = 1.6/xm, in Fig. 
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FIG. 1: Plot of the equilibrium transition temperature Te 



for aluminium as function of x — i 



Q is plotted Te(1) written in Eq. l|2*2*|l as function of 
l/l. It suggests that there is a minimum thickness be- 
low which there is no equilibrium transition tempera- 
ture. For Al, it is l^ in — 1/4036.55, and remembering 
that I = L/£ , we obtain L% in = 2.48 x 10"%, i.e. 
L min ~ 4 Angstrom 



Moreover, notice from Eq. i|19|) 
that our model is restricted to films of thickness lower 
than L% ax = 12.4£ ~ 18.4 x 10~ 6 (in Fig. Q this 
corresponds to a minimum value of x, x m i n ~ 1/12.4). 
If we take in Eq. I|14(l only the /-independent part of 
the c 4 term we would obtain u = A, i.e, no /-correction 
to the self-coupling constant. In this case, the condi- 
tion c 2 <C 1 would give a slightly higher upper limit, 
Lmax ~ 20.1£o- I n other words, taking into account the 
/-correction to the self-coupling constant constrains the 
film to be slightly thinner. Also it should be noticed that 
we have used 3-dimensional values for the coupling con- 
stants. This means that we should restrict ourselves only 
to relatively thin films; very thin films can not be phys- 
ically accommodated in the context of our model. The 
value we have obtained for £„j n (~ 4 Angstrom) cor- 
responds to 2-dimensional systems, and so, on physical 
grounds, is beyond the domain of validity of the model. 

Now, let us derive the expression of the size tempera- 
ture of the first-order transition, defined by (AT) E — 
\L(T E )/AC(T E )\, where L(T E ) is the latent heat at 
T E and AC(T E ) is the specific-heat jump. The la- 
tent heat is obtained from L(T E ) = T E AS(T E ), 
where AS(T E ) = Sqn(T e ) — Sqb(T e ) is the jump in 
the entropy (S = dT jdT) at T E between the two 
phases. Also, specific-heat jump comes from AC(T E ) = 



dSpN (Te) Sqb{Te) 



dT 



dT 



Hence, the expression for the 
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size temperature is 

(AT) J3 = T «=^°. (24) 

To obtain an numerical estimation, we can come back to 
the aluminium example and consider a sample of L w 6£o , 
i.e. near the upper limit of validity of the model. We 
get (AT)^ = 6.1 x 10~ 6 K. We see that this value is 
greater about 3.5 times than that in 1] for a bulk alu- 
minium sample. It is important to remember that (AT) 
in |]J is calculated with the specific-heat jump taken at 
the bulk second-order transition temperature To, which 
differs of a factor 1/4 from AC at the bulk first-order 
transition temperature, Tqe, as is pointed in Ref. @. 
Therefore, taking (AT) Al in calculated at Tqe yields 
(AT)q E (HLM) = (AT)q' {HLM) /A « 1.75 x 1CT 6 A, 
which jutifies the claim of an enhancement in the size 
temperature of the transition for films. 

In conclusion, we have investigated the HLM effect in 
films using a compactified version of the GL model. As a 
consequence, the free energy Ijl5(l has a logarithmic and 
linear dependence on the film thickness L through the co- 
efficients of the 4> 2 and 4 terms, respectively. However, 
it should be noticed that in spite of these differences, the 
T-dependence of Eq. (|24(l for the size temperature is for- 



mally the same as in Ref. 4J. The size temperature 
of the transition has a T - ^dependence, which means 
that the first-order phase transition is enhanced as L de- 
creases, in qualitative agreement with £|. The equilib- 
rium transition temperature of the first-order transition 
decreases with the lowering of the film thickness, imply- 
ing formally the existence of a minimal film thickness 
Lmini below which the transition is suppressed. Also, 
we notice from Eq. I|19(l that the values of L have an 
upper bound, L max = 247rA£o/e 4 - Thus the present 
model is restricted to films of thickness L such that 
Lmin < L < L max ■ For aluminium, we get respectively, 
L m l ax = 12.4& « 18.4 x 1CT 6 and L% in « 4 Angstrom. 
However, as explained in the comments below Ea. (|23ll . 
we must restrict ourselves only to relatively thin films 
which could be considered as essentially 3-dimensional 
objects. We can not on physical grounds, pretend to ap- 
ply our formalism in its present form to very thin films, 
which are too-dimensional systems. In particular, the 
value we have obtained for L m l in (Ri 4 Angstrom), is phys- 
ically meaningless, far beyond the scope of our model. 
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